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Presentation  Outline 


•  Motivation 


•  Distributions  and  distortion 


•  Examination  of  distortion  functions 


•  Numerical  example 


•  Conclusions 
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Motivation 


•  Suppose  a  risk’s  associated  severity  is  described  by  a  dis¬ 
tribution 


•  Risk  measures  summarize  distribution  (e.g.,  mean-variance 
methods,  VaR,  conditional  VaR,  distorted  expectation) 


•  Expectation  dampens  catastrophic  outcomes  -  right  tail 
may  require  further  emphasis  (risk  aversion) 


•  Questions 

-  How  do  distortions  interact  with  distributions? 

-  Which  distortion  function  and  parameters  to  select? 
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Distributions 


•  Risk  (R)  distribution 


R(x )  =  P(X  >  x 


Y  =  l)-(l-p) 


where  X 

y 

p 


severity  given  undesirable  outcome 
binary  RV  of  occurrence  (l=yes,  0=no) 
prob  of  no  undesirable  outcome 


•  Severity  (S)  distribution  (p  =  0) 


S(x)  =  P(X  >  x 


Y  =  iy 


Air  Force  Institute  of  Technology 
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Distortion  Effects 


Severity 

Figure  1.  Distortion  effects  on  exponential  distribution. 
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☆ 


Distortion 


•  Distortion  function,  g 

—  Emphasizes  worst  outcomes  (“pushes"  density  right) 

—  Forms  a  composition,  g(S(x ))  =  (goS)(x) 

—  Is  a  transformation,  g  :  [0, 1]  — ►  [0, 1] 

•  Gamma-beta  distortion  (McLeish  &l  Reesor,  2003) 

ru 

9GB (u)  =  /  -  £)6-1  exp(-t/c)dt, 

Jo 

where 

—  /  ta_1(l  -  t)b~x  exp(-t/c)dt  and  u  =  S(x ) 
j  o 

Air  Force  Institute  of  Technology  — 
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Literature  Review:  Distortion 


•  GB  family  -  six  distortions,  selected  parameters: 
gamma-beta  (GB)  (a,b,c),  beta  (c  — >  oo),  proportional 
hazard  (PH)  (6=1,  c  — »  oo),  dual  power  (DP)  (a  =  1, 
c  — >  oo),  gamma  (6  =  1),  exponential  (EX)  (a  =  1,  b  =  1) 
(McLeish  &  Reesor,  2003) 


•  Parameter  ranges:  0  <  a  <  1,  b  >  1,  and  c  >  0  are  suf¬ 
ficient  to  ensure  coherency  (i.e.,  that  risks  behave  “rea¬ 
sonably")  (Artzner,  et  al . ,  1997) 


•  Apparently  no  published  works  on  appropriate  choice  for  a 
distortion  function  or  selection  of  associated  parameters 
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Distribution  Selections 


•  Unbounded  distributions 
—  Exponential 

—  Weibull 

•  Bounded  distributions 
—  Triangular 

—  Uniform 


Air  Force  Institute  of  Technology 
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☆ 


Distortion  Function  Selections  and  Effects 


Selected  distortion  functions:  PH,  DP,  EX,  [GB] 


Table  1.  General  distortion  effects  for  survivor  function  S(x). 


Distortion 

Parameter 

{g  o  S)(x) 

Proportional  Hazard  ( gPH ) 

0  <  a  <  1 

Sa(x) 

Dual  Power  ( gDp ) 

b  >  1 

1  -  (1  -  S(x))b 

Exponential  (gsx) 

0  <  c  <  oo 

l-exp(-S(x)/c ) 
l-exp(-l/c) 

Gamma-Beta  ( gGB ) 

a,  6,  c  (as  above) 

nS(x) 

/  ta~1(l— f)6-1e_£/c  dt 

Jo 

f  fa_1  (1— t)b~1e~t/c  dt 

Jo 

Air  Force  Institute  of  Technology 
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Distortion  Parameters  for  Experimentation 


•  Background 

—  New  distortion  measure  for  GB  req’d  (expectation  N/A) 

—  If  ^  is  undistorted  median,  Rg  —  ^9°SS^^ 

•  “Region  of  sensitivity"  for  Rg  =>  1  <  Rg  <  2  (loses  track 
of  “distance  pushed") 

•  3^-factorial  design  for  GB  -  “fair"  analysis  required  each 
parameter  have  equal  influence  over  Rg  measure 

•  Face-centered  cube:  three  values,  equally  spaced 
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Selected  Distortion  Parameter  Treatments 


Table  2.  Selected  distortion  parameter  treatments. 


Distortion  (Parameter) 

Selected  Value(s) 

Rg  (%  density  shift) 

High 

0.9 

1.07  (7%) 

Proportional  Hazard  (a) 

Mid 

0.75 

1.19  (19%) 

Low 

0.6 

1.32  (32%) 

Low 

1.1 

1.07  (7%) 

Dual  Power  (6) 

Mid 

1.3 

1.19  (19%) 

High 

1.5 

1.29  (29%) 

High 

3.6 

1.07  (7%) 

Exponential  (c) 

Mid 

2.2 

1.11  (11%) 

Low 

0.8 

1.30  (30%) 

Air  Force  Institute  of  Technology 
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Analytical  Expectation  Results 


•  Explicit  expressions  for  distorted  expectation  risk  measure 


•  For  single-parameter  distortions,  numerical  results  attain¬ 
able  even  when  analytical  expectation  intractable 

Table  3.  Summary  of  risk  measures,  X  ~  exp(A). 


Distortion 

S(x ) 

E[X] 

9ph 

g— A  ax 

(A a)-1 

9dp 

1  -  (1  -  e~Xx)b 

f°°[l  —  (1  —  e~Xx)b]  dx 

JO 

9  EX 

1— exp(— e~Xx/c) 

f°°  l-exp(-e-AVc)  , 

1— exp(— 1/c) 

Jo  1— exp(— 1/c) 

Air  Force  Institute  of  Technology 
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Weibull  Distribution:  Analytical  Results 


Table  4.  Summary  of  risk  measures,  X  ~  Weib (/3,0). 


Distortion 

S(x ) 

E[X] 

gpH 

ea{-x/ey 

6  \~(  1  'I 

9DP 

1  -  (1  -  eC-*/0)'3)6 

J0°° [  1  —  (1  —  dx 

9  EX 

l-exp  (-e^/^/c) 
l-exp(-l/c) 

f°°  i— exp (— eC-^/c) 

Jo  1— exp(— 1/c) 

Air  Force  Institute  of  Technology 


Triangular  Distribution:  Analytical  Results 


Table  5.  Summary  of  risk  measures,  X  ~  tria(0i,  02, m). 


Distortion 


S(x) 


E[X] 


gpH 


(!  -  (fe-S.ff-0,))  ,  Ol<x<m 

((feT(fe-m))“’  ”i<x<02 


JTd 


(x-0i); 


(02—01  )(m— 0i) 


\  a 

J  dec 


(02— m)a+1 

(2a+l)(02-01)“ 


9  DP 


-( 

(i 


(x-0i); 


(02— 01 )  (m  — 01 ) 


(02-x)2 


(02  01 )  (02  to) 


\  b 

J  ,  6i  <  x  <  m 


\  & 

j  ,  m  <  x  <  6  2 


m  —  9\  — 


f 2  i-(i 


(m—6  i)fc+1 
(02  — 0i),J(26+l) 

(^2-x)2  \ 

(02  01 )  (02  TO.)  / 


dec 


l-exp 


(=?■ 


(g-gi)f 


c(e2-Sl)(m-0i) 


l-exp 


l-exp  (— 1/c) 
-(e2-®)2 


Oi  <  x  <  m 


rm 

Je  i 


/  -(02-x)2  \ 

\  c(e2-fii)(92-m)  / 


,  c(e2-fii)(02-m)  /  ,  ^  /i 

l-exp(-l/c)  ■  m<X<9  2 


1-exn  I  =1+ 

e  ^c(fl2-$1)(m-ei)/ 

l-exp  (-1/c) 

-(®2“x)2  ^ 


dec 


02  l-exp  ( 

Jm  1— € 


c(02-®l)(®2-m) 


— exp(— 1/c) 


dec 
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Uniform  Distribution:  Analytical  Results 


Table  6.  Summary  of  risk  measures,  X  ~  unif(0i,02). 


Distortion 

S(x ) 

E[X) 

gpH 

9dp 

(^2  0l)(6+1) 

9  ex 

l-exp(-(l-g^)/c) 

01  X  i-c+ce-^T 

1— exp(— 1/c) 

'■”2  (  l-e-i/c  J 

Air  Force  Institute  of  Technology 


Effectiveness  and  Efficiency 


•  Effectiveness’.  K  —  Rg/ \iq  =>  K  >  1 


•  Efficiency: 


K  _  change  in  fi  _  A/i 

Rg  change  in  density  Adensity 


•  Importance  of  combined  measure 

—  Without  it,  no  ability  to  distinguish  between  pairings 
with  identical  effectiveness  ( “many-to-one"  mapping) 

—  Every  increase  in  Rg  =  additional  “step"  from  SME 
recommendations  — >  undesirable  consequence 
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Effectiveness  and  Efficiency  Measures 


Table  7.  Summary  of  effectiveness  and  efficiency  measures. 


Distortion  — ► 

PH 

DP 

EX 

Measure  j 

a  =  0.9 

a  =  0.75 

a  =  0.6 

b  =  1.1 

b  =  1.3 

b  =  1.5 

c  =  3.6 

c  =  2.2 

c  =  0.8 

Exponential(3.5),  no 

=  0.285714 

0.318 

0.382 

0.476 

0.304 

0.336 

0.366 

0.306 

0.319 

0.379 

Rg 

1.072 

1.189 

1.319 

1.067 

1.188 

1.293 

1.069 

1.113 

1.303 

K 

1.111 

1.333 

1.667 

1.063 

1.177 

1.280 

1.070 

1.116 

1.327 

K/R„ 

1.037 

1.121 

1.263 

0.996 

0.991 

0.990 

1.001 

1.003 

1.019 

Weibull(2,2),  /i0  = 

1.772454 

l^g 

1.868 

2.047 

2.288 

1.845 

1.971 

2.079 

1.845 

1.891 

2.097 

Rg 

1.072 

1.189 

1.319 

1.067 

1.188 

1.293 

1.069 

1.113 

1.303 

K 

1.054 

1.155 

1.291 

1.041 

1.112 

1.173 

1.041 

1.067 

1.183 

K/R„ 

0.983 

0.971 

0.978 

0.976 

0.936 

0.907 

0.973 

0.958 

0.908 

Triangular(l,7,4), 

H  o  =  4.000 

f^g 

4.116 

4.322 

4.578 

4.103 

4.279 

4.425 

4.097 

4.159 

4.428 

Rg 

1.072 

1.189 

1.319 

1.067 

1.188 

1.293 

1.069 

1.113 

1.303 

K 

1.029 

1.080 

1.144 

1.026 

1.069 

1.106 

1.024 

1.039 

1.107 

K/R„ 

0.960 

0.909 

0.868 

0.961 

0.901 

0.856 

0.958 

0.934 

0.849 

Uniform (1,7),  /i0 

=  4.000 

l^g 

4.158 

4.429 

4.750 

4.143 

4.391 

4.600 

4.139 

4.227 

4.609 

Rg 

1.072 

1.189 

1.319 

1.067 

1.188 

1.293 

1.069 

1.113 

1.303 

K 

1.039 

1.107 

1.188 

1.036 

1.098 

1.150 

1.035 

1.057 

1.152 

K/R„ 

0.969 

0.931 

0.900 

0.971 

0.924 

0.889 

0.968 

0.949 

0.885 
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Preferred  Distortion  Functions 


Table  8.  Preferred  distortion  functions. 


Distribution 

Low  Distortion 
(0-10%) 

Moderate  Distortion 
(11-20%) 

Heavy  Distortion 
(21-30+%) 

Exponential(3.5) 

PH 

PH 

PH 

Weibull(2,2) 

PH 

PH 

PH 

Triangular(l,7,4) 

DP 

EX 

PH 

Uniform(l,7) 

DP 

EX 

PH 

•  Achieve  largest  possible  increase  in  mean  given  a  specified 
maximum  shift  in  density 


•  Shift  density  by  smallest  amount  required  to  achieve  a 
specified  increase  in  mean 
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Guidelines  for  Distortion  Selection 


•  GB:  Inability  to  analytically  compute  distorted  expecta¬ 
tion  less  appealing  choice 


•  If  SMEs  suggest  exponential,  Weibull:  PH  distortion  most 
efficient 


•  If  SMEs  suggest  triangular,  uniform:  not  as  clear 


•  If  additional  moments  are  obtained  from  distorted  dis¬ 
tribution  (e.g.  ,  cr2),  DP  and  EX  may  be  preferred  over 

PH 
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Results:  Decision  Maker  Policies 


Decision  maker’s  assigned  weights  capture  priorities 


Table  9.  Notional  point  allocations. 


Capability 

Assigned  Weight 

A 

20.0 

B 

30.0 

C 

19.0 

D 

13.0 

E 

6.0 

F 

6.0 

G 

6.0 

H 

0.0 

J 

0.0 

Air  Force  Institute  of  Technology 


Notional  Distributions 


Table  10.  Notional  distributions. 


Capability 

Distribution 

A 

Weibull(3. 5,1101) 

B 

Tria(l,  46773, 1585) 

C 

Unif(l,104) 

D 

Tria(l,104,100) 

E 

Weib(2. 04,24.73) 

F 

Weib(3. 08, 359.1) 

G 

Unif(l,100) 

H 

Exp(0.0063) 

J 

Tria(l,75,3.16) 

Air  Force  Institute  of  Technology 
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Application  of  Distortion 


•  Proposed  methodology  applies  distortion  on  distribution- 
by-distribution  basis 


•  Specific  distortions  applied  objectively  in  accordance  with 
guidelines  previously  discussed 


Table  11.  Selected  distortion  function  application  results. 


Capability 

Distribution 

(1  -p)no 

Rg 

Combination 

(1  ~P)Hg 

A 

Weib(3. 5,1101) 

24.766 

1.20 

PH,  a  =  0.735 

27.043 

B 

Tria(l, 46773, 1585) 

32.239 

1.30 

PH,  a  =  0.62 

42.641 

c 

Unif(l,104) 

37.504 

1.19 

EX,  c  =  1.3 

42.264 

D 

Tria(l,104,100) 

33.670 

1.13 

EX,  c  =  1.9 

37.189 

E 

Wei  b(2. 04, 24. 73) 

7.887 

1.06 

PH,  a  =  0.915 

8.238 

F 

Weib(3. 08, 359.1) 

9.631 

1.06 

PH,  a  =  0.915 

9.913 

G 

Unif(l,  100) 

18.938 

1.06 

DP,  b  =  1.09 

19.737 

H 

Exp(0.0063) 

38.095 

1.0 

N/A 

38.095 

J 

Tria(l,75,3. 16) 

13.193 

1.0 

N/A 

13.193 

Air  Force  Institute  of  Technology 
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☆ 


Linear  Programming  Formulation 


9  6 

Maximize  S{ 

2=1  Jl  =  l 
6 

subject  to  kj  x j  =  25 

j=l 

0  <  ay  <  1,  j  =  1,2, 

where 

is  risk  expectation  accompanying  Capability  i 
rrii  j  denotes  mitigation  to  Capability  i  by  system  j 
Xj  is  “amount”  of  each  mitigator  to  be  purchased 
kj  is  cost  of  any  complete  mitigator  j  (25  unit  budget) 

Air  Force  Institute  of  Technology 
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Optimal  Purchase  Plan  (LP  Solution) 


Table  12.  Notional  acquisition  recommendations. 


Risk  Measure 

Mitigator 

1 

2 

3 

4 

5 

6 

None 

1.0 

1.0 

0.2 

0.0 

0.0 

1.0 

Undistorted  Expectation 

1.0 

1.0 

0.0 

0.0 

0.25 

1.0 

Distorted  Expectation 

1.0 

1.0 

0.0 

0.25 

0.0 

1.0 

Air  Force  Institute  of  Technology 
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Areas  for  Further  Study 


•  Computing  expectation  of  multi-parameter  distortions 


•  Measures  other  than  Rg  should  be  considered 

•  Study  of  correlation  between  Pearson’s  skewness  coeffi¬ 
cient  and  normalized  mean 


•  Application  of  distortion  functions  to  other  distributions 


•  Effects  of  distortion  on  variance 

Air  Force  Institute  of  Technology 
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Open 


•  Comments 


•  Questions 


Forum 


Air  Force  Institute  of  Technology 
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Coherency 


Artzner,  et  al.  (1997) 

i.  Sub-additivity. 

p(x  +  Y)  <  p(X)  +  p(Y) 

ii.  Homogeneity. 

p(t-X)  =  t-  p(X) 

iii.  Monotonicity. 

p(X)>p(Y),  \f  X  <Y 

iv.  Risk-free  condition. 

p(X  +  r  •  n)  =  p(X)  —  n 

Air  Force  Institute  of  Technology 


Graphical  Effects:  Weibull 


Severity 

Figure  3.  Relative  freq  density  for  severity,  Weib(2,2)  distribution,  given 
distortion  parameters  a  =  0.6,  b  =  1.5,  c  =  0.8  (solid  is  no 

distortion, - GB,  •  •  ••  PH,  -  • DP,  —  •  —  EX). 
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Typical  n  Effects  Plot  (Weibull) 


Coded  Distortion  Parameter  Value 


Figure  4.  Expected  value  versus  coded  distortion  parameters,  Weib(2,2) 
distribution,  given  distortion  parameter  ranges 
a  =  [0.525,0.975],  b  =  [1, 1.6],  c  =  [0.1,  4.3]. 


Air  Force  Institute  of  Technology 


Graphical  Effects:  Triangular 


Severity 


Figure  5.  Relative  freq  density  for  severity,  tria(l,7,4)  distribution,  given 
distortion  parameters  a  =  0.6,  b  =  1.5,  c  =  0.8  (solid  is  no 

distortion, - GB,  •  •  ••  PH,  -• —  DP,  —  •  —  EX). 
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Figure  7.  Relative  freq  density  for  severity,  unif(l,7)  distribution,  given 
distortion  parameters  a  =  0.6,  b  =  1.5,  c  =  0.8  (solid  is  no 

distortion, - GB,  •  •  ••  PH,  —  •  —  DP,  —  •  —  EX). 
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